Comment on "Probing vortex unbinding via dipole fluctuations" 



O 

o 

(N 
< 



43 

o 



-i— > 



I 

G 

O 

o 



(N 
> 

m 

o 
o 

-I— > 



I 

c 

o 
o 



x 



S. E. Korshunov 

L. D. Landau Institute for Theoretical Physics, Kosygina 2, Moscow 119334, Russia 

(Dated: December 5, 2006) 

We demonstrate that the method suggested by Fertig and Straley [Phys. Rev. B 66, 201402 
(2002)] for the identification of different phases in two-dimensional XY models does not allow one 
to make any unambiguous conclusions and make a tentative proposal of another approach to this 
problem. 
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Recently, Fertig and Straley [l[ have proposed a new 
method for the identification of different phases in two- 
dimensional XY models. They have introduced the so- 
called extended dipole moment [2| -P e " t (where a = x,y), 
which generalizes the notion of vortex gas dipole moment 
for the systems with periodic boundary conditions, and 
have suggested studying the probability of large fluctu- 
ations of this quantity. According to Ref. [l|, T{L), the 
probability of finding P e ^ t (or, equivalently, P^ xt ) equal 
to (n + 1/2)L (where L is the linear size of the system 
and n is an arbitrary integer) , must have different depen- 
dences on L (when L — > oo) in phases with bound and 
unbound vortices. 

With the present Comment, we would like to point 
out that the asymptotic behavior of J-(L) has to be the 
same \T{L) oc 1/L] independent of whether the vortices 
are bound in pairs or free, and therefore the analysis of 
J-{L) does not allow one to distinguish different phases. 

Consider a phase with unbound vortices. In the first 
approximation, one can assume that in such a phase 
the vortex positions are completely uncorrelated. In that 
case, it is evident from the definition [l], 0] of P " t 
that this quantity must have uniform distribution, and 
therefore T(L) = 1/L . The interaction between vor- 
tices makes the integer values of P" xt /L more preferable 
than half-integer, and therefore leads to some decrease of 
T(L). ^ 

Consider now a phase in which strong interaction be- 
tween vortices binds them into small neutral pairs well 
separated from each other. When typical distance be- 
tween vortex dipoles is much larger than their size, their 
interaction is weak, which suggests that one can ne- 
glect the correlations between the orientations of differ- 
ent dipoles. In such an approximation, the total dipole 
moment of the system, P a = YljPf' ^ s the sum of N 
independent random variables p" (N being the number 
of vortex pairs), so for N 3> 1 and \P\ <§; N it has Gaus- 
sian distribution with the width a oc TV 1 / 2 . Therefore, 
the probability to have P a = ±L /2 is given by 
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P a = P c " t (mod L), Eq. (p} also gives the dependence of 
?{L) on L. 

Now we have to return to the interaction between vor- 
tex dipoles. If it would be dependent directly on the rel- 
ative orientation of the dipoles, it would lead to suppres- 
sion (or enhancement) of the fluctuations of total dipole 
moment. However, in two dimensions the dipole-dipolc 
interaction 
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where we have taken into account that in large systems 
N = cL 2 , c being the concentration of vortex pairs. Since 



depends not on the relative orientation of two dipoles Q 
but on their orientations with respect to the line which 
connects them. This interaction does not force the to- 
tal dipole moment of a diluted gas of dipoles to become 
smaller (or larger), and therefore cannot lead to suppres- 
sion (or enhancement) of its fluctuations in comparison 
with the case of nonintcracting dipoles. 

In Ref. [D, the fluctuations of dipole moment in the 
phase with bound vortices were analyzed with the help 
of the duality transformation with a subsequent replace- 
ment of a discreet solid-on-solid (SOS) model by a con- 
tinuous one (with a renormalized coupling). In terms of 
the original XY model, this corresponds to replacing the 
cosine interaction by a harmonic one. Naturally, such a 
replacement leads to the complete suppression of vortices, 
which fixes the total dipole moment at zero [as follows 
also from Eq. (8) of Ref. [l[ . It is clear that the disregard 
of vortices makes this approach utterly unsuitable for an- 
alyzing the fluctuations of vortex gas dipole moment. 

Since we have found that the asymptotic dependence 
of F(L) on L is the same both for noninteracting vortices 
and for strongly bound vortices, one can expect than it 
will be the same also in all intermediate cases. There- 
fore, the numerical calculation of T(L) cannot help one 
distinguish the phases with bound and unbound vortices. 

The observation in the numerical simulations of Ref. 1 
of the situations in which T[V) decays with the increase 
of L much faster than 1 / L can be explained as a manifes- 
tation of transitional regimes in which the main contribu- 
tion to T{V) is related to the formation of large dipoles 
and not to the reorientation of small dipoles. However, 
our analysis suggests that with a further increase of L, 
a crossover to the regime in which the mechanism de- 
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scribed in this Comment is the dominant one must occur. 
Note that in the case of strong coupling, the coefficient 
between T[V) and 1/L is exponentially small in vortex 
pair concentration c, and therefore the observation of the 
asymptotic dependence of T(V) requires rather large val- 
ues of L (with \xiL ^> 1/c). 

Possibly, a more consistent approach to the identifi- 
cation of different phases in two-dimensional XY mod- 
els can be based on analyzing in dynamic simulations 
the behavior of "world lines" of vortices in the three- 
dimensional space-time. These world lines are in many 
respects analogous to vortex lines in a three-dimensional 
XY model. In particular, they have to be continuous 
due to the conservation of the topological charge. It fol- 
lows from the results of Ref. that in terms of vortex 
world lines, the main difference between the phases with 
bound and unbound vortices consists in the absence (or 
existence) of the world lines of infinite length (crossing 
a "whole sample" in time direction). In the phase with 
bound vortices all world lines have to form closed loops, 
whereas in the phase with unbound vortices the concen- 
tration of infinite world lines (which can be associated 
with free vortices) has to be nonvanishing. 

This suggests that in the framework of a dynamic 



description, the concentration of free vortices is a well 
defined quantity (in contrast to thermodynamic simula- 
tions, in which it is impossible to introduce an algorithm 
for counting the number of free vortices in a given vortex 
configuration). Therefore, a numerical calculation of this 
concentration may turn out to be a useful instrument for 
distinguishing the phases with and without free vortices. 

The author is grateful to H.A. Fertig for useful ex- 
change of opinions and to V.V. Lebedev for enlightening 
discussion. 
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